Consequences of statistical sense determination for WIMP directional detection 
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We study the consequences of limited recoil sense reconstruction on the number of events required 
to reject isotropy and detect a WIMP signal using a directional detector. For a constant probability 
of determining the sense correctly, 3-d read-out and zero background, we find that as the probability 
is decreased from 1.0 to 0.75 the number of events required to reject isotropy using the mean angle 
statistic is increased by a factor of a few. As the probability is decreased further the number of 
events required using this statistic increases sharply, and in fact isotropy can be rejected more 
easily by discarding the sense information and using axial statistics. This however requires an 
order of magnitude more events than vectorial data with perfect sense determination. We also 
consider energy dependent probabilities of correctly measuring the sense, 2-d read-out and non-zero 
background. Our main conclusion regarding the sense determination is that correctly determining 
the sense of the abundant, but less anisotropic, low energy recoils is most important for minimising 
the number of events required. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 



Weakly Interacting Massive Particle (WIMP) di- 
rect detection experiments aim to directly detect non- 
baryonic dark matter via the elastic scattering of WIMPs 
on detector nuclei [l| , and are presently reaching the sen- 
sitivity required to detect the lightest neutralino (which 
in most supersymmetry models is the lightest supersym- 
metric particle and an excellent WIMP candidate) . Since 
the expected event rates are very small ( 0(1O~ 5 — 1) 
counts kg _1 day _1 ) distinguishing a putative WIMP sig- 
nal from backgrounds due to, for instance, neutrons from 
cosmic-ray induced muons or natural radioactivity, is 
crucial. The direction dependence of the WIMP scat- 
tering rate (due to the Earth's motion with respect to 
the Galactic rest frame) [2( provides a potential WIMP 
'smoking gun' and low pressure gas time projection 
chambers (TPCs), such as DRIFT (Directional Recoil 
Identification From Tracks) 0,11] and NEWAGE [J], see 
also Ref . [|| , seem to offer the best prospects for a detec- 
tor capable of measuring the directions of sub-100 keV 
nuclear recoils. 

Early studies found that in principle as few as 30 events 
would be required to distinguish a WIMP induced signal 
from isotropic backgrounds [7|, |8| . In reality the number 
of events, and hence the exposure required, depends on 
the detector properties including the energy threshold, 
background event rate, whether the read-out measures 
the recoil momentum in 2 or 3 dimensions (and if 2-d 
in which plane) [t| E3, EH, 03 and whether the sense 



(i.e. the absolute sign + f> or — p) of the recoil momentum 
vectors can be measured [1, Eo, E2, El, EH . 

The factor with the biggest effect on the number of 
events required is whether or not the sense of the re- 
coils can be measured d, E3, EH • If the data is axial (i.e. 
the sense can not be measured) the exposure is increased, 
compared to the vectorial case, by one order of magnitude 
for 3-d read-out and at least two orders of magnitude for 
2-d read-out in the optimal plane [23[ [l2j]- At low ener- 
gies, the energy deposition dE/dx of nuclear recoils, and 
hence the density of ionisation created, is predicted to 
slowly decrease with decreasing recoil energy 14 1. Thus 
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the sense of a recoil track is, in principle, measurable by 
determining the direction in which the ionisation density 
decreases along the track. 

Previous work on simulating recoil tracks using 
SRIM2003 15] indicated that the ionization density dis- 
tributions reconstructed from recoil tracks were close to 
uniform due to both fluctuations in the production of ion- 
isation and diffusion during the drift of this ionisation to 
the read-out plane [1, EH ■ It was therefore not possible to 
determine the absolute signs of the reconstructed recoil 
vectors in these simulations. It should be noted however 
that while SRIM2003 predicts sulfur recoil ranges and 
quenching factors (fraction of recoil energy going into 
ionisation) in agreement with experimental data [161 ], it 
was not designed to model recoils in gaseous targets. 

Whether or not the sense of low energy nuclear re- 
coils can be determined is therefore an important ques- 
tion which needs to be resolved experimentally. Ongo- 
ing studies [13, EH indicate that the expected decrease 
in dE/dx is observable, but fluctuations in the ioniza- 
tion distribution caused by the ionization process and 
diffusion may wash out the information in the ioniza- 
tion distribution leading to errors in sense reconstruc- 
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tion. In other words, if the sense of a recoil of en- 
ergy E is determined via a hypothesis test on some 
parameter M derived from its reconstructed ionization 
distribution, then the PDFs <7(M|forward sense, E) and 
g(M | backward sense, E) will not be cleanly separated, i.e. 
the significance level, S(E), is greater than zero and the 
power, V(E), is less than one. Here, the significance 
level is the probability of rejecting the forward sense hy- 
pothesis when it is true, and the power is the probabil- 
ity of rejecting the backward sense hypothesis when it is 
false [lj|. Provided S(E) < 0.5 and V{E) > 0.5, the cor- 
rect recoil sense will be reconstructed more often than not 
with probability P(correct sense) E) = 1 — S(E). In this 
Brief Report we therefore study the effect of this proba- 
bilistic sense reconstruction on the number of events re- 
quired to distinguish a WIMP induced recoil signal from 
isotropic backgrounds. 



II. CALCULATIONS 

We use the same statistical techniques and methods 
for calculating the directional nuclear recoil spectrum as 
in Refs. [§, [l(| G3. We briefly summarise these pro- 
cedures here, for further details see these references and 
Ref. [2(| ■ We consider the simplest possible model for the 
Milky Way halo, an isotropic sphere with local density 
p = 0.3 GeVcm -3 and a Maxwellian velocity distribution 
with dispersion <r v = 270 km s _1 , and fix the WIMP mass 
at m x — 100 keV. Our simulated detector is a TPC filled 
with 0.05 bar CS 2 200 fim pitch micropixcl readout 

plane, a 10 cm drift length over which a uniform electric 
field of 1 kV cm -1 is applied, and is based on the design of 
the DRIFT-I/II detector Q. We use the SRIM2003 [H 
package to generate sulfur recoil tracks and for 3-d read- 
out recoil directions are reconstructed as the principal 
axis of the charge distributions recorded by the pixels. 
2-d read-out would measure the projection of the recoil 
momentum vector into a plane fixed on the Earth. The 
degree of anisotropy of the 2-d recoil angles (and hence 
the detectability of a WIMP signal) depends on the ori- 
entation of the read-out plane [lOj, [ll|, [l2j . Here we focus 
on the optimal case of a read-out plane with normal per- 
pendicular to the spin axis of the Earth. No simulations 
of the an gula r resolution function of a 2-d detector are 
available [24 ] we therefore assume perfect resolution in 
this case (and hence the resulting numbers of events are 
lower limits on the number which would be required in 
reality) . For primary recoil energies below E t h = 20 keV 
the short track length (3-4 pixels) and multiple scattering 
make it impossible to reconstruct the track direction in 
our simulated detector, we therefore only consider events 
with energies above this threshold. Zero background is 
the goal of the next generation of experiments made from 
low activity materials with efficient gamma rejection and 
shielding, located deep underground [2l|, however we 
investigate the effect of non-zero isotropic background 
by varying the ratio of the background and signal event 



rates. 

Recoil directions in 3-d and 2-d constitute points on 
the unit sphere and circle respectively. For 3-d data the 
most powerful test for rejecting isotropy uses the average 
of the cosine of the angle between the direction of solar 
motion and the recoil direction, (cos 9) [9( . For 2-d data 
the most powerful test is the Rayleigh test [22| which 
uses the mean resultant length of the projected recoil 
vectors which, modulo fluctuations, should be zero for 
data drawn from an isotropic distribution [ic| . 

We calculate the probability distribution of the rele- 
vant statistic, for a given number of events N, by Monte 
Carlo generating 10 5 experiments each observing N re- 
coils drawn from our simulated 3-d or 2-d distributions. 
We then compare this with the null distribution of the 
statistic, under the assumption of isotropy and calculate 
the rejection and acceptance factors, R and A, at each 
value T of the statistic. The rejection factor is the prob- 
ability of measuring a value of the statistic less than T 
if the null (isotropic) hypothesis is true or equivalently 
the confidence with which the null hypothesis can be re- 
jected given that measured value of the statistic. The 
acceptance is the probability of measuring a value of the 
statistic larger than T if the alternative hypothesis is true 
or equivalently the fraction of experiments in which the 
alternative hypothesis is true that measure a larger ab- 
solute value of the test statistic and hence reject the null 
hypothesis at confidence level R. Clearly a high value 
of R is required to reject the null hypothesis, while a 
high A is also required, otherwise any one experiment 
might not be able to reject the null hypothesis at the 
given R or the null hypothesis might be erroneously re- 
jected. We therefore find the number of events required 
for A C = R C = 0.95, N 95 . 

We consider several functional forms for the energy de- 
pendence of the probability of correctly determining the 
recoil sense, P(correct sense|P) (hereafter P CS (E)). Our 
simplest model assumes an energy independent proba- 
bility, P CS (E) = PloOkcVi With 0.55 < PlOOkeV 

< 1.0. 

In reality it is likely that it will be easier to determine 
the sense of higher ener gy recoils due to their longer 
track lengths (c.f. Ref. [3). We therefore also con- 
sider linearly increasing probability P CS (E) = aE + b. 
We parametrise this function in terms of the values of 
Pes (Eth = 20keV), where P t h is the threshold energy as 
above, and P CS (E = 100 keV) = PiookeV- We consider 
two parameter sets for this parametrization, P cs (£'th) = 
0.50 with 0.55 < PiookcV < 1-0, and P cs (E th ) = 0.75 
with 0.75 < PiookcV < 1-0, the latter being more op- 
timistic about reconstruction at low energies. In each 
case for E > Eu m , where P cs (Pii m ) = 1.0, we set 
P CS (E > E Hm ) = 1.0. 



III. RESULTS 

For 3-d read-out, zero background and constant sense 
determination probability the number of events required 
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FIG. 1: The dependence of the number of WIMP events 
required to reject isotropy at 95% confidence in 95% of ex- 
periments, N95, using the value of {cos 8) on the probabil- 
ity of correctly determining the recoil sense. The circles are 
for energy independent probability, P CS (E) = PiookcV, tri- 
angles (squares) for the probability linearly increasing from 
P ca (E th = 20keV) = 0.75 (0.5) to P CS (E = 100 keV) = PiookcV 
The dotted line shows the number of events required using the 
value of (| cos#|), which does not require the sense of the re- 
coils. The panels are (top row from left to right and then 
bottom row) zero background and signal to (isotropic) back- 
ground ratio S/N — 10, 1 and 0.1. In all cases, including 
the axial {|cos#|) statistic, we assume 3-d readout with an 
energy threshold of E t h = 20keV and take into account the 
uncertainty in the reconstruction of the recoil direction. 
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FIG. 2: As fig.[TJfor 2-d read-out using the Rayleigh statistic. 
As discussed in the text the angular resolution is not taken 
into account in this case, and these numbers are hence lower 
limits on the number of events required by a real detector. 



to reject isotropy using the (cos 9) statistic initially de- 
creases fairly slowly as the probability is decreased (from 
N 95 = 11 for FiookoV = 1.0 to N 9S = 39 for FiookoV = 
0.8). The increase becomes more rapid as PiookcV is de- 
creased further and N95 — > 00 as PiookeV — * 0.5. The 
linearly increasing sense determination probability with 



P cs (E th ) = 0.75 and PiookeV = 1.0(0.75) requires ~ 3 
times more (the same number of) events as the con- 
stant probability with the same value of PiookeV- For 
■Pes (-Eth) = 0.50 the increase is more dramatic (more 
than an order of magnitude for PiookcV = 1-0). In 
fact in this case (and also for the constant probability, 
if PiookcV <^ 0.7) fewer events are needed to reject 
isotropy using the axial (| cos0|) statistic which does not 
use any information about the sense of the recoil. This 
is at first glance surprising, however in these cases the 
low values of P cs , especially at low energies where the re- 
coil rate is higher and the anisotropy in directions lower, 
significantly decrease the anisotropy in the observed re- 
coil momentum vectors and hence the sensitivity of the 
(cos 6*) statistic. This does not affect the (|cos0|) statis- 
tic, as it is independent of the value of P cs , being sensitive 
to the concentration of recoil axes, rather than vectors, 
around the line of solar motion. The number of events 
required using the axial statistic is, however, an order 
of magnitude larger than for vectorial data with perfect 
sense determination. As the signal to noise ratio is de- 
creased the numbers of WIMP events required in each 
case are increased, by factors of 1.1 — 1.2, 2 — 4 and ~ 10 
for S/N — 10, 1 and 0.1 respectively. 

For 2-d read-out with perfect angular resolution, the 
numbers of events required are a factor of roughly 2 larger 
than for the same sense determination probability and 
signal to noise ratio for 3-d read-out with the recoil direc- 
tion reconstruction uncertainty taken into account. We 
caution that projection effects will make angular resolu- 
tion a more significant factor for 2-d read-out than for 
3-d read-out. Estimates based on the projected length 
of recoil tracks in the planes indicate that the required 
number of events would increase by at least a further 
factor of 2 [Io|. 

Finally we examine whether it is more important to 
correctly determine the sense of the abundant low en- 
ergy recoils or the rare, but more anisotropic, high en- 
ergy recoils. To do this we consider step function sense 
determination probabilities: P CS {E < E stcp ) — 0.5(1.0), 
P CS {E > E stcp ) = 1-0(0.5) for 3-d read-out and zero 
background. These correspond to the sense being un- 
determined (perfectly determined) for low energy recoils 
and perfectly determined (undetermined) for high energy 
recoils. It should be noted that neither of these possibil- 
ities, in particular the later, are particularly physically 
plausible. As the energy above which the sense can be 
(perfectly) determined, E siep , is increased the number of 
events required increases rapidly (by an order of magni- 
tude as .Estep is increased from 20keV to 60keV). Con- 
versely for the case where the sense of high energy recoils 
can not be determined as -E s tep is decreased Ng$ increases 
slowly initially, and then more rapidly for E siep < 40 keV. 
This indicates that determining the sense of the common, 
but less anisotropic, low energy events is most important 
for minimising the number of event required to reject 
isotropy. 
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IV. SUMMARY 

In this Brief Report we have investigated the effect of 
probabilistic recoil sense determination on the number of 
events required by a directional detection experiment to 
reject isotropy, and detect a WIMP signal. We have con- 
sidered a constant sense determination probability and 
also increasing probability with increasing recoil energy 
(higher energy recoils have longer tracks and hence it 
is likely that it will be easier to determine their sense). 
For a constant probability of determining the sense cor- 
rectly, 3-d read-out and zero background, we find that 
as the probability is decreased from 1.0 to 0.75 the num- 
ber of events required to reject isotropy using the mean 
angle statistic is increased by a factor of a few. As the 
probability is decreased further the number of events re- 
quired using this statistic increases sharply, and in fact 
for PiookeV < 0.7 isotropy can be rejected more easily by 
discarding the sense information and using axial statistics 
(this does however require an order of magnitude more 
events than vectorial data with perfect sense determina- 
tion) . 

The linearly increasing sense determination probabil- 
ity with P cs (E th ) = 0.75(0.50) and Pi 00kc v = 1.0 re- 



quires ~ 3 times (an order of magnitude) more events 
than the constant probability with the same value of 
-PiookeV- This suggests that maximising the probabil- 
ity of correctly determining the sense of the low energy, 
more common, but less anisotropic events is most impor- 
tant for minimising the number of events required. We 
have confirmed this conclusion by considering step func- 
tion sense determination probabilities. 

We also considered non-zero background and 2-d read- 
out. As the signal to noise ratio is decreased the numbers 
of WIMP events required in each case are increased, by 
factors of - 1.1 - 1.2, 2 - 4 and 10 for S/N = 10, 1 and 
0.1 respectively. For 2-d read-out with perfect angular 
resolution, the numbers of events required are roughly a 
factor of ~ 2 larger than for 3-d read-out with the recoil 
direction reconstruction uncertainty taken into account 
(with the same sense determination probability and sig- 
nal to noise ratio). 
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